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Abstract 

In this paper we consider a connection between the famous Skorohod embedding 
problem and the Shiryaev inverse problem for the first hitting time distribution 
of a Brownian motion: given a probability distribution, F, find a boundary such 
that the first hitting time distribution is F. By randomizing the initial state of 
the process we show that the inverse problem becomes analytically tractable. 
The randomization of the initial state allows us to significantly extend the class 
of target distributions in the case of a linear boundary and moreover allows us 
to establish connection with the Skorohod embedding problem. 
Key words: First hitting time, Shiryaev's inverse boundary problem, 
Skorohod embedding 

1. Introduction 

Let Xt be an arbitrary process with the left-continuous sample paths and 
let b{t) be a continuous absorbing boundary satisfying the condition Xq > 6(0). 
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The random variable 

inf{t >0:Xt< b{t)}, if there exists t such that Xt < b{t), 
oo, otherwise, 
is caUed the first hitting time for the process Xt. The problem 

Given a process, Xt, and the boundary, b(t), (1) 
find the distribution Fr{t) = P (t < i). 

is the starting point of a very rich research area in the theory of stochastic 
processes. If Xt is a diffusion process, the problem of finding the distribution of 
T is a classical one. The first papers on the problem were published by P. Levy, 
A. Khintchine and A. Kolmogorov in the 1920s. These results are discussed 
in the monograpfp] |Khintchine (1933). Khintchine gave a complete solution of 



the problem for sufficiently smooth boundaries. This solution was expressed 
in terms of a boundary value problem for the associated partial differential 
operator (infinitesimal generator). Since then many books and research papers 



were published in this area. The monograph Lerche ( 1986 1 summarizes known 



analytical results obtained by the mid-1980s. In Karatzas and Shreve (2005) the 
link between analytical methods and the martingale approach is considered, the 
paper 



Durbin ( 1971 ) discusses computational aspects of the problem, while the 



Taylor expansions of the probability distribution of r are considered in |Hobson| 



et al. (1999). 



The following, inverse to ([!]), problem was proposed by Albert Shiryaev in 



^Khintchine called |TJ the second problem of diffusion. 
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the mid-TO's in his Banach Center lectures: 



Given a process, Xt, and a distribution, F{t), (2) 
find a boundary, b{t), such that V {t < t) = F(t). 

Problem ([2]), to the best of our knowledge, was proposed by A. Shiryaev in the 
case when Xt is a Brownian motion and F{t) is an exponential distribution. If Xt 
is a Brownian motion with a random initial value, i.e., the process Xt — S, + Wt 
is a Brownian motion starting at a random point ^, we will refer to Problem ([2| 
as a generalized Shiryaev's problem (GSP). 

The third problem considered in this paper is the Skorohod embedding prob- 
lem: 

Given a probability measure, fi, with a finite second moment, (3) 
and a Wiener process, Wt, find an integrable stopping time, r*, 
such that the distribution ofWr, is /i. 

This problem has a long list of references. Here we mention only the original 



paper Skorohod (1965) and a survey by Obloj (2004) where many important 
applications and various solutions of the Skorohod embedding problem can be 
found. 

It is intuitively appealing that there must be a connection between the prob- 
lems ([2]) and Indeed if a monotone continuous function solves Prob- 
lem ([2| for a Wiener process Wt, and the distribution function F{t), then the 
first hitting time, r*, is a solution of the Skorohod problem with P {Wr, ^ t) = 
F (b^^{t)y The only problem with the latter condition is that is usually 
unknown in analytical form. Moreover, the boundary might not belong to the 
class of monotone functions. 

Models using the first hitting time distribution find extensive application in 
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the areas of portfolio credit risk modeling (see Iscoe and Kreinin (19991, Iscoe 



et al. (1999), Schmidt and Novikov (2008)) and pricing of credit derivatives 



(Avellaneda and Zhu (2001), Hull and White (2001)). In these contexts the 



process represents the so-called distance to default of an obligor (see Avellaneda 



and Zhu (2001)), while the first hitting time represents a default event. The 
boundary therefore acts a barrier separating the healthy states of the obligor 
from the default state. For this reason, the boundary, b{t), is often called the 
default boundary in the applied literature on credit risk modeling. In particular, 
an interesting model of default events with a randomized boundary was proposed 
Schmidt and Novikov] ( |2008[ ). 



The inverse problem (|2| was considered in Iscoe et al. (1999) when the 
process Xf is a Brownian random wall([^ A detailed analysis of the inverse 



problem in the discrete time setting is given in Iscoe and Kreinin (1999) as 
well as a Monte Carlo based solution. Their approach is applicable to a much 
more general class of processes Xt, not just Brownian random walks, and is 
computationally simple to implement. 

Existence of the solution to the continuous-time inverse problem ^ is an- 



alyzed in Chadam et al. (2006). In Peskir (20021, an integral equation for the 



boundary is derived when Xt is a Brownian motion. A general analysis of the 
integral equations for the boundary and existence and uniqueness theorems are 



considered in Jaimungal et al. ( 2009a I . The randomized inverse problem is con- 



sidered in Jaimungal et al. ( 2009b I and in the preliminary publication Jackson 



et al. (2009). 



The Shiryaev problem is notoriously difhcult and analytical solutions are 



Shepp 


( 


1967 


), 


Salminen (1988 


1, 


Lerche 



•^A Brownian random walk is a discrete time process with Gaussian increments and variance 
proportional to the time step. 
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(1986), Peskir and Shiryaev (2006), Breiman (1967), Alili and Patie (2005), 



Novikov (1981), Peskir (2002)). However, existence of the solution to the prob- 



lem has been proven for an arbitrary target distribution, F, by [Dudley and| 



Gutmann 


(1977 


) and by 


Anulova ( 


1980 



tence theorem in its most general form cannot be proven for the the Generalized 



Shiryaev's problem and in fact counter-examples do exist (see Jaimungal et al. 



(2009b) and Proposition 2 below). Nonetheless, the randomized version does 
admit closed form solutions for a large class of distributions of the first hit- 



ting time, F, including the gamma distribution (see Jackson et al. (2009) and 



more generally Jaimungal et al. (2009b)) as well as a subset of one-sided stable 
distributions (see Corollary 1 below). 

One of our main goals is to establish a connection between the Generalized 
Shiryaev's problem and the Skorohod embedding problem. The remainder of 
this paper is organized as follows. In Section[2]we proved a new short derivation 
for the Laplace transform of the distribution of the first hitting time and the 
distribution of the initial random position ^ of the process Xt . When the bound- 
ary is linear, we find solutions for a class of mixtures of gamma distributions 
of the first hitting time and for a class of stable distributions. In Section |3] we 
analyze the structure of the solutions to the Shiryaev's problem and introduce 
the minimal solution which possess a very elegant structure. The randomization 
of the initial state of the process allows us to stretch the boundary and trans- 
form it into a straight line. In this case we have a simple relation between the 
distributions of r and Xr- This observation allows us to connect the Skorohod 
embedding problem to the Generalized Shiryaev problem as shown in Section [4j 

This paper is self-contained; it represents an extended version of the talk 
given at the 5*'* Bachellier colloquium, January 2011, Metabief. 
Acknowledgement. We are very grateful to Yuri Kabanov and Tom Salisbury 
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for pointing out the appealing similarity of the Problems under consideration 
and asking about the connection between the inverse Problem ([2| and the Skoro- 
hod embedding problem. Wc would also like to thank Lane Hughston and Mark 
Davis for the interesting comments on earlier stages of this work. The authors 
are indebted to Albert Shiryaev, Alexander Novikov, Raphael Douady, Isaac 
Sonin and Nizar Touzi for the interesting discussions during the 5*'' Bachellier 
colloquium. 

2. Randomization of the initial state 

2.1. General equation 

Let {n,^,dt,P) be a filtered probability space and let 17 = C ([0,oo)) x 
be the Cartesian product of the space of continuous functions on the positive 
semi-axis and M^. Consider a random variable ^ on this probability space. We 
assume that {^t} = {S^J^} V 23(]R^) is an augmented filtration where is the 
natural filtration, generated by the standard Wiener process Wt and assume 

that *8 = [J dt- We also assume that the random variable ^ is independent of 

t>o 

the process Wt- 

Let Xt — ^ + Wt and let P (0 < ^ < +oo) = 1. Consider now the linear 
boundary, b{t) — kt, k > Q. Denote the first hitting time of the process Xt to 
b{t) byrW, i.e., 

T^'') = inf{i : i > 0, C + Wt<kt}. 

The Generalized Shiryaev's Problem (GSP) for the process Xt is formulated as 
follows: 

Given a distribution F{t), find ^, 

(4) 

such that P (r^'') < = F{t). 

The triplet (^, t'-'^-' , k) is called a solution of the Generalized Shiryaev's Problem. 
We are interested in the solutions to the GSP and their properties. In 
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particular, we are interested in the solution of Problem Q when F{t) belongs 
to the class of Gamma distributions, F G ^x.-y, with the probability density 
function 

p^{t) = A • ^ ; • exp(-At), 7 > 0, A > 0. 



Remark 1. The case 7 = 1 corresponds to the original formulation of Shiryaev's 
problem with a randomized initial point. In this case the distribution F{t) = 
1 - e-^\ t > 0. 



To solve the problem, we will derive a connection between the Laplace trans- 
forms of the target distribution F(t) and the initial starting point ^. To this 
end, denote 



f{s) = E[e— ] 



dF{t), 



s e 



(5) 



where, as usual, M-)_ — {s : s > 0}, is the set of non-negative real numbers, and 
let g{s) = E[e^^^. The function /(s) on [0,oo), is completely monotone (see 



Feller 



(1971)), f em, where 



(F-V(s) 

9Jt=<iy(s):(-l)"^^>0, n = 0,l,..., seM+ 



The class of completely monotone functions form an algebra: the sum and 
the product of completely monotone functions belong to 931. According to the 
classical Bernstein's Theorem (see 



Feller 



(1971)), if /(s) e at and /(O) = 1, 



then there exists a cumulative distribution function, F{t), satisfying ([s]). 

The following statement for completely monotone functions is very well 



known (see Feller (1971), Criterion 2) 



Proposition 1. // G{s) is a completely monotone function, u{s) > and the 
first derivative u' G 971 then G{u{s)) G 9Jl. 

Let us now derive the main equation for the Laplace transforms of the distribu- 
tions of T and ^. 
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Theorem 1. If a random variable ^ is a solution to the Problem Q then the 
function g{s) satisfies the equation 

9{s) = f{sk + s^/2). (6) 

Proof. The process Xt is a martingale with respect to the fihration {5't}f>o- 
Consider the exponential martingale Mt = exp (— sXj — s'^t/2), (s > 0). It is 
not difficult to show that if ^ is finite almost surely then the exponential mar- 
tingale is bounded for all t < t. Therefore, by the optional stopping theorem, 

E[A/^] = E[Mo]. 

Moreover, the expected value 

E[Mo] =E[e-^«]. 
Furthermore, he have Xr = kr and 

Mr = exp (-r • [sk + s'^ /2)) . 

Therefore E[Af^] = f{ks + s'^/2). Finally, we obtain Equation ■ 

2.2. Gamma-distributed first hitting time 

In the case / e ^x,-^, the random variable ^ solving Problem Q admits the 
following simple probabilistic interpretation. 

Theorem 2. Let f{t) ~ p-y{t), for some 7 > 0. Suppose that 2 A < k^ . Then 
g{s) satisfies 

_ (fc + VP-2A)T (fc - VP - 2A)T 

^ ~ {s + k + Vfc2 - 2A)'' ' {s + k~Vk^- 2Xy ■ 

and the random variable 

e = 6+6, (8) 

where ^1 and ^2 o,fe independent gamma-distributed random variables with a 
common shape parameter 7; ^1 ^ r(7, Ai) and ^2 ^ ^(7, A2), with 

Ai = fc - Vfc^ - 2A, A2 = fc + Vk^ - 2A. (9) 

Proof. We have 

A'T 



From (|6| we obtain 



2A 



s2 + 2ks + 2A 
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If fc^ > 2A the quadratic equation 

+ 2ks + 2X = 

has real roots — Ai and — A2. Then, taking into account that A1A2 = 2A, we find 

'^'^ = (a^) • (a^) ■ 

The additive representation ([8| for ^ follows immediately from the latter equa- 
tion. I 

Remark 2. It follows from Theorem[2]that if the GSP has a solution for k = k^, 
given f(t) e Ta^^, then it also has a solution for any k > k^,. We will generalize 
this property for a larger class of distributions in the next section. 

Theorem [2] can also be generalized in the following direction. Consider a 
class of random variables Mr which are mixtures of gamma-distributed random 
variables with respect to shape 7 and scale A parameters. Write, 

T = Mix(rA,^,Ai) , 

where ^ is a mixing measure having a support, (A, 7) ^ S — [0, A*] x [0, 00), 
/X [S) = 1, and A* = P /2. Then 

/»:=E[e-]=/^(^dM(A,7). 



In this case 



s 



where Ai(A) and A2(A) satisfy ([9]). Therefore, we see that the random variable 
^ which leads to the hitting time distribution Mix (r>,,-y, /i) is in fact a mixture 
of the sums of independent gamma-distributed random variables. 

This generalization is interesting because it is known that the class of finite 
mixtures of gamma-distributed random variables 77 G Mr form a dense set in 
the space of random variables. Moreover, an arbitrary random variable 77* can 



be obtained as a weak limit as the number of gamma distributions n — >■ oo, 

Vn — >V*: J7„ e Mr. 
Unfortunately, this sequence of approximations can not be used to extend the 



solution ( 10 1 through weak convergence because the parameter A has a bounded 
support. In particular, we cannot approximate constant almost surely random 
variables using the random variables, ?]„ G Mr- 
2.3. Stable distribution of the first hitting time 

So far, we have considered the case when r has finite mean, Et < +cxi. 
However, it is also interesting to analyze the case of one-sided stable distribu- 
tions with infinite expected value. For this case, we take the slope coefficient 
A; = 0. Then, if one takes the random variable r such that 

E[e-"^] = e-=^°, oO, 0<a<l, 

then it follows from (|6| that 

g{s) = exp (-cis'^) , 
where (3 = 2a and ci — i^. Thus, we obtain 

Corollary 1. Ifr has a stable distribution with the parameter a, < a < 1/2, 
then ^ has a stable distribution with the parameter 2a. 

3. Structure of the solutions to Shiryaev's problem 

3.1. General results 

Theorem [2] describes the solutions to the Shiryaev's problem when the first 
hitting time t is gamma-distributed. In this case, it can be shown that the 
condition k > k^, — \/2 A is necessary and sufficient for existence of the random 



variable ^ (see Jackson et al. (20091 and Jaimungal et al. (2009b)). Thus if the 
problem has a solution for the slope fc* then it has a solution for any fc > fc*. 
The following statement generalizes this result. 
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Theorem 3. Given t, suppose there exist > and — ^(fc*) such that 
Equation Q is satisfied: 

g4s) :=E[e-^«-] = f {hs + 3^2) . 

Then for any k > k^, there exists a random variable S,{k) such that 

E [e-'^«('^-)] = / (fcs + s V2) , s > 0. (11) 

Proof. The function g{s) = / (A:s + satisfies the relation g(0) = 1. There- 

fore, it is enough to prove that g £ DJl. Let k > k^,. Consider the function 



Obviously, u 



„2 



u{s) = v/s2 + 2ks + kl - fc*. 
and 



— +ks = — ^ h k^,u{s), s > 0. 



The first derivative of u(s) is 
d 



ds 



(s) = u'{s) = {s + k) ■ {s^ + 2ks + kl) 



Let us prove that the function u'{s) e 9Jt. We have 

■u'(s) • (m(s) + fc*) = s + /c. 
Then if fc > /c* , we derive from ( |12[ ) the inequality 

u'{s) > 1. 

Denote n*'' derivative of the function u{s) by m(")(s). (n = 1, 2, . . . ). Differen- 
tiating Equation ( 12 ) we find 

. {u{s) + k,) + {u'{s)f = 1. (14) 



(12) 



(13) 



The latter Equation and ( 13 1 imply the inequality 

u^'^\s)<0, s>0. 



Denote m — [?^/2]. From Equation (14) we derive 

m 



(15) 
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The coefficients Cn,j satisfy tfie relations 
Cn+1,1 = n = 2,3, 

Cn+i,j = Cn,j + Cn,j^ij n = 2, 3, . . . , j = 2, . . . , m — 1. 
Taking into account tliat C2,i = 1, we find 

C„,, = (^^j, j-l,2,...,m-l. (16) 

The coefficient C„ m satisfies the relation 



r if ^ is odd, 

' ' * otherwise, 



2 Vm/ ' 



Thus, Cn,j > in Equation (15). Then from (15) we prove by induction that 

(-l)"+iu(")(s) >0 for all n = 1,2,.... 
Therefore u'{s) E SJl. Then we have 

Since S 9Jt and u'{s) e 9Jt we derive from Proposition [I] that g{s) e 9Jl, as 
was to be proved. B 
The next proposition demonstrates that there exist parameters of the GSP 
such that the problem does not have a solution. 

Proposition 2. Suppose E[t^] < oo. Then the Problem Q) does not have a 
solution if k < ■\/E[r]/(T(T). 

Proof. Indeed, from Equation ^ we find 

dgjs) ^ df{ks + sy2) 
ds ds 

and 



{k + s), 



d^gjs) _ dV(fc,s + ,sV2) , d/(fc,s + ,sV2) 

d^^ ^^^'> + d^ 

Substituting s = into these equations, we find 

E[^] = k-E[T], 
E[^2] = /fc2E[r2] -E[r]. 
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The equations for the first two moments of imply 

k^a^ir) ^ E[t] + C7^{^) > E[t] 

Finally, we obtain 



air) ■ 

Proposition [2] is thus proved. B 

Remark 3. Propositioiij2] provides only a necessary condition for existence of 
a solution to Problem In particular, if r ~ r(7. A), then the necessary 
and sufficient condition for existence is fc > -s/ZA and the minimal solution 
corresponds to fc* = %/2 A. 

Remark 4. If r is a constant with probability 1, Equation ^ does not have a 
solution. The GSP can not be solved using our randomization approach in this 
case. 

Theorem [3] and Proposition [2] imply that the admissible set for the coefficient 
k is either the semi-infinite interval [fc*,cxD) or an empty set. If t has finite 
first two moments and the random variable ^ exists for some fc* > then for 
each fc > fc* there exists a random variable ^(fc) solving Problem Q. In this 
manner, we can obtain a family of solutions S = {C(^)} parameterized by the 
slope coefficient fc. Furthermore, the variance of the random variable ^(fc) is 
a monotone function of fc and the one with minimal variance is what we call 
the minimal solution to Problem Q. In particular, if r has an exponential 
distribution with parameter A, then the random variable ^* corresponds to the 
minimal solution has slope fc* = \/2A. In this case the random variable ^(fc*) is 
Erlang distributed with order two. 

3.2. Esscher families 

Consider a non-negative random variable 77 and denote fri{s) = E[e^'"'], s > 
its Laplace transform. It is convenient to introduce a family of random vari- 
ables, 

£(7?) = {7^ : E[e-^""] - /,(a + s)/Ua), a > o}. 
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which we call the Esscher family generated by the random variable rj. As before, 
let r be the first hitting time and /(s) — E[e^*'^] be its Laplace transform. 
Suppose the generalized Shiryaev's Problem Q has a solution ^ = ^(fco) for 
some fixed slope fco- Consider now two Esscher families £(t) and £(^) generated 
by T and ^, respectively. 

Proposition 3. For each k > ko there exist a = a{k) and j3 — f3{k) such 
that Ta G £(t) and € £(^) and {Ta,£,p,k) is a solution to the Generalized 
Shiryaev's problem. 

Proof. Denote A = fc - fco, (A > 0), and take a = A • fco + /3 = A. 

Consider the random variable € £(0- Then after simple transformations we 
find 

. , , /(a + (fco + A)-g + .V2) f{a + ks + sy2) . 

9p{s) = J— = J—. = /a(fcs + S ). 

Thus, the triplet (tq, f^, fc) form a solution of the problem. I 

4. Connection to Skorohod problem. 

We are now in a position to construct a solution to the Skorohod problem 
through the solution to the Generalized Shiryaev's problem. Suppose that the 
distribution of the stopped process has a non-negative support, P(Xt- > 0) = 
1, and the continuous distribution, Fx{x) = V{Xr < x), can be matched by a 
corresponding distribution of the random variable ^ for some fc > 0. Taking 
into account that Xt- = kr for the linear boundary b{t) = kt, (fc > 0), we 
immediately obtain that P(r < t) = Fx{kt) and the Laplace transform of the 
density of r satisfies 

In this case 
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Similarly, the Skorohod problem can be solved for the distributions with V{Xt- < 
0) = 1. In this case, the first hitting time is understood as 

T = infji : Xt> kt, k <0}, 

i.e., the first time the process touches the boundary fi'om below. Clearly, the 
initial point ^ must then have support on the negative real axis. 

The solution to the general Skorohod problem - without restrictions on the 
support of the distribution of X^- - can then be obtained as follows: Let us 
represent the random variable X^ as a mixture of the random variables X^ > 
and X_ < 0: 

{X+, with probability p+=¥ (X^ > 0), 
X_ , with probability p_ = 1 — p+ . 

Suppose that the Skorohod problem for the random variables X^ and X_ can 
be solved by the boundaries 6+ (t) — k+ ■ t and 5_ (t) — k^ ■ t and the random 
variables ^-|_ and respectively (see Figure [l]). Then we have the following 
Proposition 4. The random variable 

^ with probability p^, 

1 with probability 

and the boundary b{t) = 6+ U 6_ solve the Skorohod problem for the random 
variable X^.. 

Remark 5. In this manner, the Skorohod problem is solved by randomizing the 
starting point of the Brownian motion over the entire real line, and searching 
for the first hitting time when the process enters the wedge region described by 
the lines fe+ t and k-t. 

5. Conclusion 

We have demonstrated that randomization of the initial state of the process 
is a very powerful tool for solving and tying together the Generalized Shiryaev's 
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g+(^) 



ai 



Figure 1: Randomizaton in the Skorc 



bJt)=lct 



anct^tjjc Skorohod Embedding Problem. The randomizat ion of the 
in^iajl sljate of the proceSs^aJlows us to reduce the problem to the linear bound- 
iand obtain closed- form sof(!»tiQnJ^fe(^verJ^itportant cla4ses of the 

n of the firatWpiJrti^iijtLe. 
lineariB(ii..ti£)rif of the boundaVSfvna,kes the reldMon b^T^We^ the Bhiryaev's 



d 8ko]'ohod problems transparent. Furtnifimilfc, it allows us to develop a new 



to the Skordhod probleiti. 



We close this paper with a brief discussion of the directions for future re- 
search. In the present paper we discussed the case of a scalar process Wt whose 
distribution of the first hitting time must be matched. One interesting open 
question is to prove existence of the solution to the Generalized Shirayaev's 
problem when r has a continuous infinitely divisible distribution. Another in- 
teresting generalization that requires much attention is the Generalized Shi- 
rayaev's problem for a vector- valued process Wt- Such generalization will prove 
extremely useful in application settings as well as providing an interesting math- 
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ematical playground. 
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